Sieving hydrogen based on its high compressibility 
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Abstract 

A molecular sieve for hydrogen is presented based on a carbon nanotube intramolecular junction and 
a Cgo- The small interspace formed between C^o and junction provides a size changeable channel for 
the permselectivity of hydrogen while blocking Ne and At. The sieving mechanism is due to the high 
compressibility of hydrogen. 
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With growing demand on clean energy, the production and separation of hydrogen with great 
amount represent a topic of intense interests. One of the straightforward and commercial ways 
to separate hydrogen is to use a specifically designed molecular sieve with high selectivity and 
permeability. Efforts on using the traditional molecular sieve, such as polymeric membrane, have 
been ongoing for some time and the progress has been made along this direction. yj,|2|] However, 
the major problem related to polymeric membranes is the strong trade-off between selectivity and 
permeability, yj, 120 Thus, the search of new sieves or new methods to improve the permselectiv- 
ity of hydrogen from mixtures is much needed to match the growing demand in the hydrogen 
energy. [I Ij 

In molecular sieves, the major mechanism for the filtration usually bases on the difference in 
size, shape or polarity of varies molecules in the mixture. H 2 has almost spherical shape and no 
polarity, thus the probable sieving mechanism for H 2 is limited to use its small size. Attempts to 
design sieves based on the size difference have been made over the past few years.O, |4f| This kind 
of sieve is quite useful for molecules with large size difference. However, for molecules having 
comparable sizes with H 2 , the efficient way for the separation is still lack of. 

The purpose of this letter is to find an alternative mechanism to design sieves for hydrogen. 
We notice that hydrogen has a novel property, namely the highly compressible character, which 
is absent for most other molecules, o, |6y A hydrogen molecule has only two electrons on fairly 
compact Is orbitals, the van der Waals repulsion between two H 2 molecules is much weaker than 
typical, so that it is highly compressible all the way up to and above GPa pressure. The soft 
repulsion is also found for H 2 interacting with other systems, such as graphite. [29] The highly 
compressible character makes hydrogen have much more chances to diffuse at confinement space 
over other gases. Furthermore, H 2 is the lightest molecule in nature, at a given temperature, so 
H 2 could have higher mobility in confinement space than other gases. Considering the above two 
features, it is possible to design a molecular sieve to filtrate hydrogen. 

The ideal candidate for current studies could be carbon nanotubes (CNTs) due to its tubular 
structures and small pore sizes. In fact, over the past few years, many studies have focused on the 
investigation of the selectivity and transport properties of light gases in CNTs. [7 - 111 However the 
use of CNTs as a molecular sieve for hydrogen is still a great challenge. If one hopes to directly 
use the pore of CNTs, the CNT should have comparable pore size with H 2 (~ 3A). Such small pore 
size almost approaches to both theoretical predicated limitation and the smallest CNT fabricated 

Up to now, as far as we know, the reliable design based on CNTs is still 



in experiments. 



lack of. 

To avoid the above mentioned limitation, we do not directly use the channel of CNTs, but the 
interspace formed between C^q and CNTs. The inset of Fig. 1(a) presents the geometry of the 
molecular sieve. It is made of an intramolecular junction and a Cm. To form the intramolecular 



junction between two tubes, the method presented in Ref. 11151- 1911 is adopted. Combining the 



well-deve 



oped experimental techniques for producing j unctions ||20l- l22|l and for inserting C^q into 



CNTs|23|-|25|], our design could be easily fabricated. In current studies, the left tube is (20,0) CNT, 
which is larger than C 6 o, while the right one is (5,5) CNT, smaller than C 60 . The interspace formed 
between C 60 and CNTs could provide a nano channel for the separation of hydrogen using its high 
compressibility. 

To compare the permselectivity of hydrogen with other gases, two typical gases are also consid- 
ered, which are Ne and Ar. All the molecules/atoms have almost spherical shape and no polarity, 
the major difference is their size. Among these molecules, Ar is the one with largest diameter 
(about 3.8A); while H 2 and Ne have the similar sizes of about 3.3A and 3.1 A respectively. Be- 
cause the difference in sizes is not prominent and H 2 even has larger size than Ne, the separation 
could be very difficult, even impossible in traditional molecular sieves. If hydrogen can be sepa- 
rated from those molecules mentioned above, this design could also work for most other types of 
molecules, such as CH 4 , N 2 and 2 . It needs to be pointed out that, all the three molecules can 
easily pass both (20,0) and (5,5) CNTs. 

To validate this molecular sieve, several calculations have been carried out. First, to check the 
size and stability of the narrow space formed between C^o and CNT, we have studied the absorp- 
tion behavior of C^o- Second, to verify if a single H 2 molecule can easily pass the sieve without 
the highly compressibility emerged, the diffusion barrier for one molecule passing the sieve has 
been calculated. Third, to validate the importance of the high compressibility, a comparing study 
on the transport properties for each pure gas in the sieve is carried out. Finally, the permselectivity 
of H 2 mixed with other gases is also calculated. This work demonstrates that, if the highly com- 
pressible character of hydrogen are taken into account, our design does work as a molecular sieve 
for hydrogen. 

We have addressed the above questions using molecular dynamic (MD) simulation. To cal- 
culate the permselectivity and flux, the molecular sieve with length about 36 A is dipped into a 
reservoir, as shown in Fig. 1(b). Two plates, labeled as A and B, are placed at the two boundaries 
of the reservoir perpendicular to the molecular sieve. These plates, made of the corresponding gas 



molecules, are rigid with face-centered cubic (111) planes. In all simulations, plate B is fixed, 
while plate A is movable. In all studies, the system temperature keeps at 300 K by Nose-Hoover 



thermostat. (12611 The equations of motion are integrated by using the predictor-corrector algorithm 
with time step of 1 fs. Periodic boundary conditions are only applied in directions perpendicular 
to the CNTs. The pressure difference around C^ is defined as AP = P le f t - Pnghu where Pj e f, is 
fluidic pressure of the reservoir, and P r i g ht keeps zero in all simulations. In current studies, AP is 
changed by placing the plate A at different positions. For each move of plate A, the large CNT is 
blocked by fixing a few atoms/molecules inside it for first 0.4 ns so that the fluid in the reservoir is 
fully equilibrated. Afterwards, all gas atoms/molecules are movable, and the system is simulated 
for about 0.6 ns to determine the flux of gases. 

The interatomic potential between carbon atoms is the Tersoff-type many-body potential with 



parameters given by Brenner. [27 



281] The weak interaction between carbon atoms in CNTs and 
C(,o is modeled by the Lennard- Jones (L-J) potential, V vc i W = Ci 2 /r 12 -C 6 /r 6 with C6=20 eVA 6 and 
C\i - 2.48 xlO 4 eVA 12 . Interactions between H 2 molecules are modeled by the Silvera-Goldman 
potential, which correctly reproduce the equation of states for H 2 in wide range of pressures. [5Q 
The weak interaction between H 2 and carbon is described by the recently fitted potential by some 
of US. I29I which reproduces the results over a wide range of repulsive and attractive regions as 



calculated by high level ab-initio methods. [30] The interactions of Ne and Ar are L-J potential with 
the parameters taken from reference[31]. Within the pressure up to 1 GPa, the L-J potential could 
accurately predicate melting curve and other properties which are comparable to experimental 



results .1131 



3211 The interactions of C-Ne and C-Ar are also L-J potential and calculated using 



Lorentz-Berthlot combination rules with parameters of C-C taken from Ref.[33H. 

Without filling any gases, there is one stable adsorption site for C^q around the intramolecular 
junction. Fig. 1(a) shows the potential energy of the system as C(,o moving along axial direction. 
At this stable adsorption site, C^o is right located at the axis center near the junction. The minimal 
distance between and the junction is around 3.3 A as indicated in the inset of Fig. 1 (a) 
labeled as 'interspace A'. The interlayer distance from C(,o to tube wall is about 4.3 A marked as 
'interspace B'. The stable adsorption position is a reflection of the good geometric match between 
junction and Ca), namely, favorable interactions between pairs of C atoms on CNTs and C(,q reach 
maximum. The adsorption energy is about 1.17eV lower than that for Ceo far from the junction. 
The stable adsorption guarantees C(, staying near the junction, which is important for the current 
sieve. Due to the weak interaction between C 60 and tube wall, C 6 o can be moved around the stable 



position in a certain range, which forms a 'size changeable channel'. 

The diffusion barriers for single molecule/atom passing the molecular sieve are calculated, 
which is plotted in Fig. 2(a). We find that, when the molecule locates at interspace A, the global 
potential maximum appears. As the molecule locates at interspace B (right above Ceo), there is 
another local potential energy maximum (flat), which is more obvious for Ar. There is a strong 
correlation between the position of molecule/atom and displacement of the mass center of C 60 . 
Fig. 2 (b) and (c) show the associated motion of the mass center of Ceo as the molecule moving 
along the tube axial direction. As the molecule approaches to the interspace B, C 60 is pushed 
backward and begins to move to the tube wall. The displacement of C 60 produces increasing 
space for the molecule. It reaches the maximum as the gas molecule arrives interspace A. At 
the same time, the potential energy reaches the maximum. After the molecule passes interspace 
A, Ceo is ready to move back to its equilibrium position. This implies that the energy barrier 
is mainly due to the displacement of Ceo- The larger the displacement of Ceo is, the higher the 
energy barrier is. For Ar, the displacement of C 60 along radial/axial direction from the equilibrium 
position is around 1.61A/1.71A, while for H 2 or Ne, the corresponding value is 1.35A/1.07A and 
1.36A/1.08A, respectively. The diffusion barrier for H 2 and Ne has the similar value around 0.44 
eV (0.45eV for Ne and 0.43eV for H 2 ), which is lower than the barrier of 0.68eV for Ar. The 
diffusion barrier suggests that, for a single molecule/atom, it is very difficult to pass the molecular 
sieve at room temperature. As we will discuss below, if the high compressibility of hydrogen is 
emerged (plenty of H 2 molecules appear), the difference in transport behavior among H 2 , Ne and 
Ar will be significant. 

In Fig. 3(a), we have shown the flux via AP for pure H 2 , Ne and Ar at 300 K separately. 
From this figure, one can see that, with the increase of pressure difference up to 1 GPa, only H 2 
molecules are excluded out from the molecular sieve, while other gases (Ne and Ar) are blocked 
by Ceo with zero permeation. As we have discussed that Ne has almost the same diffusion barrier 
as H 2 (see figure 2(a)). Thus, H 2 and Ne should have comparable chance to pass the sieve if the 
diffusion of the single atom/molecule dominates this process. However the permeability of Ne is 
much lower than that of H 2 . The mystery can be disclosed in terms of the high compressibility of 
hydrogen. It is well known that the diffusion rate is proportional to the inverse of the square root of 
molecular mass. However the mass issue is not enough to explain the transport difference between 
H 2 and Ne. On the contrary, the high compressibility of H 2 plays a key role. Comparing with Ne, 
this character is not quite obvious when few molecules appear or the pressure is low (figure 2 (a)). 
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With the increase of pressure, the high compressibility (both soft repulsion of H 2 - H 2 and H 2 - C) 
becomes more and more important. 

The high compressibility makes H 2 have more chance to accumulate in interspace B. We have 
calculated the average number of molecules located within a narrow window around interspace B. 
This narrow window has the width of 4 A along the axial direction and centered on C(, . Fig. 3 
(b) shows the average number of molecules/atoms within this narrow window via AP. One can 
see that, the number of H 2 in this area is much more than Ne. And the number increases with the 
increase of AP. For Ne, even at 1 GPa, the number is less than two, which is several times smaller 
than that of H 2 . No Ar atom is found at interspace B even at 1 GPa. The molecules located at the 
interspace B will not only increase the chance moving to interspace A, but also lower the diffusion 
barrier due to the inward motion of C^ driven by the molecules located around interspace B. With 
growing number of H 2 moving into interspace B, H 2 is more easily to overcome the energy barrier 
and be excluded out. 

To further verify that transport properties of hydrogen is dominated by the high compressibility, 
rather than the small mass of H 2 , we have carried out two additional calculations, which is shown 
in the inset of Fig. 3 (a). In the first run (long dashed line), the mass of H 2 is artificially enlarged 
ten times. In this case, we found that the flux is simply reduced according to the change of mass. 
In the second run (dot-dashed line), the interaction of H 2 -H 2 and H 2 -C is replaced by a typical L-J 
potential. J34I1 Our calculations show that, in this case, the transport of hydrogen is terminated. 

The permselectivity of H 2 mixed with Ne and Ar is further studied separately. In this calcula- 
tion, the number ratio is chosen as H 2 : M (M=Ne or Ar)=l:l, 3:1, 9:1. Fig. 4(a) and (b) shows 
the flux of two components mixed gases as a function of AP. One can see that, in all studied cases, 
the permeation of H 2 is nonzero and increase with the increase of AP. In all the simulations, there 
is no Ar atom found to pass the molecular sieve. In the mixture of H 2 and Ne, the flux of Ne, as 
the pressure difference raised up to around 0.6 GPa, becomes nonzero. However it is still much 
smaller than that for H 2 (Please see the Supporting Information for a animation of this process). 
Detailed studies show that, the leakage of a few Ne molecules is mainly 'abducted' by H 2 . That 
is, by chance, some Ne are surrounded by H 2 , these Ne atoms pass the molecular sieve associ- 
ated with H 2 . For the mixture of three types of gas molecules {H 2 , Ne and Ar) with the ratio of 
1:1:1, our calculations also show that, only H 2 pass through the molecular sieve when AP is in the 
range 0.2GPa ~ 0.6GPa, as shown in Fig. 4 (c) (Please also see the Supporting Information for 
a animation of this process). Similarly a few of Ne also leak out from the molecular sieve as AP 
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raising up to 0.6 GPa. As AP up to lGPa, no Ar penetrates the molecular sieve during at least Ins 
simulation time. 

The final question, which should be addressed, relates to the requirements on using the high 
compressible character of H 2 . For this purpose, we have checked the sieving function for hydrogen 
in a few combinations of intramolecular junction and C 60 . Our results show that the position of 
Ceo and the size of the interspace formed between the tube wall and Ceo are the important factors. 
First, the size of interspace should not be too small, otherwise the leakage of H 2 is blocked. The 
size of interspace should also not be too large, otherwise other gases could penetrate and the high 
compressibility of H 2 plays less role. We found that the interspace between the tube wall and C 6 o 
should be around 3.5A, correspondingly, the radius of the larger CNT in the sieve around 16A is 
the ideal one. For the smaller tube, it needs to have enough space for hydrogen passing. Second, 
the potential well of Ceo should not be too stiff, otherwise Ceo could not move to produce more 
space. For example, the intramolecular junction formed between (20,0) and (10,0), which fills the 
requirements, also has the function as a molecular sieve. We believe that, besides of the current 
design, there could be many structures for using the high compressible character and the small 
mass of H 2 . 

In conclusion, by taking the high compressibility of hydrogen into account, we have proposed a 
potential molecular sieve for hydrogen based on C 60 and carbon nanotube intramolecular junction. 
Our calculations show that the narrow interspace formed between the tube wall and C 60 provides 
a size changeable channel. The diffusion barriers for a single molecule/atom passing this channel 
are quite high (0.45eV for Ne, 0.43eV for H 2 , and 0.68eV for Ar). The diffusion barriers imply 
that, for a single molecule/atom, it is difficult to pass the molecular sieve at room temperature. If 
plenty of H 2 appear ( the high compressibility of H 2 is emerged), the transport behavior of H 2 , 
Ne and Ar is significant difference. Using molecular dynamics simulations, we have studied the 
permselectivity of H 2 /Ar, H 2 /Ne and H 2 /Ar/Ne mixtures. Our calculations show that no Ar is 
found to pass the sieve, while H 2 can pass the sieve with noticeable flux. A few Ne atoms can 
also leak out, however the amount is much smaller than H 2 . The selective permeation of H 2 is the 
result of the high compressibility of H 2 . The current design could be an alternative technique for 
the separating of H 2 from most of mixtures. 
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FIG. 1. (Color online) (a) The potential energy change of the molecular sieve as C(,o approaching to the 
junction. The initial axial distance of C(,o (the original of x-axis) is around 17A far away from the junction. 
The energy has been shifted so that the energy at initial position of Ceo is zero. Inset: The top and horizontal 
views of the molecular sieve, which is composed with a carbon nanotube intramolecular junction and an 
encapsulated Ceo- (b) The illustrations of the simulation system for calculating the flux/permeability of 
gases, where the plate B is fixed, while plate A is movable to tune the pressure difference around the sieve. 
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FIG. 2. (Color online) The potential energy change of system (a), the radial (b) and axial displacement (c) 
of C6o as the gas molecule moves along axial direction. The initial axial distance of molecule (the original 
of x-axis) is around 20A far away from the junction. The energy has been shifted so that the energy at initial 
position of molecule is zero. 
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FIG. 3. (Color online) (a): The average flux varied as the function of fluid pressure difference in the single 
component gases of H2, Ne and At. Inset: The flux for hydrogen in two different cases. In the first case 
(long dashed line), the mass of H2 is artificially enlarged ten times. In the second (dot-dashed line), the 
interaction of H2-H2 and H2-C is replaced by a typical L-J potential. For comparing, the flux for hydrogen 
in normal condition is also shown (solid line), (b): The average number of molecules located in interspace 
B. 
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FIG. 4. (Color online) The average flux varied as the function of fluid pressure difference in the multi- 
component gas molecules. (a) Mixed fy/Ne with number ratio of 1:1/3:1/9:1. (b)Mixed Hx/Ar with number 
ratio of 1:1/3:1/9:1. (c)Mixed H 2 /Ne/Ar with number ratio of 1:1:1. 
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